Search for: All records

Abstract Dust-obscured galaxies (DOGs) with extremely red optical-to-infrared colors are often associated with intense starburst and active galactic nucleus (AGN) activity. Studying DOGs can provide insights into the processes that drive the growth of galaxies and their central supermassive black holes. However, the general DOG population is heterogeneous, spanning a wide range of evolutionary stages, and has X-ray obscuring column densities (N_H) covering low to high levels. In this work, we focus on seven high Eddington ratio DOGs ( $\log {\lambda }_{\mathrm{Edd}}\gtrsim -0.5$) to examine their X-ray obscuration properties using new and archival X-ray observations. We confirm that these systems are generally heavily obscured, with six out of seven havingN_H ≳ 10²³cm⁻²and three out of seven havingN_H ≳ 10²⁴cm⁻². Based on the observed similarity with the rare Hot DOG population, we argue that both high-λ_EddDOGs and Hot DOGs likely trace the postmerger phase, during which AGNs are enshrouded by large columns of dust-rich material. 

Abstract Binary supermassive black holes (SMBHs) are consequences of galaxy mergers and dominate the low-frequency gravitational-wave background. Finding binary SMBHs in existing time-domain observations has proven difficult, as their periodic, electromagnetic signals can be confused with the natural variability of single quasars. In this work, we investigate the effects of host-galaxy contamination and survey design (cadence and duration) on the detectability of binary SMBHs with the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We simulate millions of LSST light curves of single and binary quasars, with a distribution of quasar and host-galaxy properties motivated by empirical observations and the anticipated LSST detection space. We then apply simple sinusoidal curve fits as a potential computationally inexpensive detection method. We find that host-galaxy contamination will increase false-positive rates and decrease binary parameter recovery rates. Lower-mass, lower-luminosity binary systems are most likely to be negatively affected by host-galaxy contamination. We also find that monitoring duration affects binary detection more than survey effective cadence for this detection method. As the light-curve duration increases, false-positive rates are suppressed and binary parameter recovery rates, especially for binary periods, are improved. Increasing the light-curve duration from 5 to 10 yr shows the most dramatic improvement for successful binary detection and false-positive rejection, with additional improvement from extending the light-curve duration to 20 yr. The observation duration increase is especially critical for recovering binary periods that are longer than a decade. 

Abstract We present a systematic investigation of extremely X-ray variable active galactic nuclei (AGNs) in the ≈5.3 deg²XMM-SERVS XMM-LSS region. Eight variable AGNs are identified with rest-frame 2 keV flux density variability amplitudes around 6–12. We comprehensively analyze the X-ray and multiwavelength data to probe the origin of their extreme X-ray variability. It is found that their extreme X-ray variability can be ascribed to changing accretion state or changing obscuration from dust-free absorbers. For five AGNs, their X-ray variability is attributed to changing accretion state, supported by contemporaneous multiwavelength variability and the absence of X-ray absorption in the low-state spectra. With new Multiple Mirror Telescope (MMT) spectra for four of these sources, we confirm one changing-look AGN. One MMT AGN lacks multiepoch spectroscopic observations, while the other two AGNs do not exhibit changing-look behavior, likely because the MMT observations did not capture their high states. The X-ray variability of the other three AGNs is explained by changing obscuration, and they show only mild long-term optical/IR variability. The absorbers of these sources are likely clumpy accretion-disk winds, with variable column densities and covering factors along the lines of sight. 

Abstract We present photometric selection of type 1 quasars in the ≈5.3 deg²XMM-Large Scale Structure survey field with machine learning. We constructed our training and blind-test samples using spectroscopically identified Sloan Digital Sky Survey quasars, galaxies, and stars. We utilized the XGBoost machine learning method to select a total of 1591 quasars. We assessed the classification performance based on the blind-test sample, and the outcome was favorable, demonstrating high reliability (≈99.9%) and good completeness (≈87.5%). We used XGBoost to estimate photometric redshifts of our selected quasars. The estimated photometric redshifts span a range from 0.41 to 3.75. The outlier fraction of these photometric redshift estimates is ≈17%, and the normalized median absolute deviation (σ_NMAD) is ≈0.07. To study the quasar disk–corona connection, we constructed a subsample of 1016 quasars with Hyper Suprime-Cami < 22.5 after excluding radio-loud and potentially X-ray-absorbed quasars. The relation between the optical-to-X-ray power-law slope parameter (α_OX) and the 2500 Å monochromatic luminosity (L_2500Å) for this subsample is ${\alpha }_{\mathrm{OX}}=(-0.156\pm 0.007)\log L_{2500\AA }+(3.175\pm 0.211)$with a dispersion of 0.159. We found this correlation in good agreement with the correlations in previous studies. We explored several factors, which may bias theα_OX–L_2500Årelation, and found that their effects are not significant. We discussed possible evolution of theα_OX–L_2500Årelation with respect toL_2500Åor redshift. 

Abstract The coevolution of supermassive black holes and their host galaxies represents a fundamental question in astrophysics. One approach to investigating this question involves comparing the star formation rates (SFRs) of active galactic nuclei (AGNs) with those of typical star-forming galaxies. At relatively low redshifts (z≲ 1), radio AGNs manifest diminished SFRs, indicating suppressed star formation, but their behavior at higher redshifts is unclear. To examine this, we leveraged galaxy and radio-AGN data from the well-characterized W-CDF-S, ELAIS-S1, and XMM-LSS fields. We established two mass-complete reference star-forming galaxy samples and two radio-AGN samples, consisting of 1763 and 6766 radio AGNs, the former being higher in purity and the latter more complete. We subsequently computed star-forming fractions (f_SF; the fraction of star-forming galaxies to all galaxies) for galaxies and radio-AGN host galaxies and conducted a robust comparison between them up toz≈ 3. We found that the tendency for radio AGNs to reside in massive galaxies primarily accounts for their lowf_SF, which also shows a strong negative dependence uponM_⋆and a strong positive evolution withz. To investigate further the star formation characteristics of those star-forming radio AGNs, we constructed the star-forming main sequence (MS) and investigated the behavior of the position of AGNs relative to the MS atz≈ 0–3. Our results reveal that radio AGNs display lower SFRs than star-forming galaxies in the low-zand high-M_⋆regime and, conversely, exhibit comparable or higher SFRs than MS star-forming galaxies at higher redshifts or lowerM_⋆. 

Abstract Supermassive black holes (SMBHs) can grow through both accretion and mergers. It is still unclear how SMBHs evolve under these two channels from high redshifts to the SMBH population we observe in the local Universe. Observations can directly constrain the accretion channel but cannot effectively constrain mergers yet, while cosmological simulations provide galaxy merger information but can hardly return accretion properties consistent with observations. In this work, we combine the observed accretion channel and the simulated merger channel, taking advantage of observations and cosmological simulations, to depict a realistic evolution pattern of the SMBH population. With this methodology, we can derive the scaling relation between the black hole mass (M_BH) and host-galaxy stellar mass (M_⋆), and the local black hole mass function (BHMF). Our scaling relation is lower than those based on dynamically measuredM_BH, supporting the claim that dynamically measured SMBH samples may be biased. We show that the scaling relation has little redshift evolution. The BHMF steadily increases fromz= 4 toz= 1 and remains largely unchanged fromz= 1 toz= 0. The overall SMBH growth is generally dominated by the accretion channel, with possible exceptions at high mass (M_BH≳ 10⁸M_⊙orM_⋆≳ 10¹¹M_⊙) and low redshift (z≲ 1). We also predict that around 25% of the total SMBH mass budget in the local Universe may be locked within long-lived, wandering SMBHs, and the wandering mass fraction and wandering SMBH counts increase withM_⋆. 

Abstract Tidal disruption events (TDEs) could be an important growth channel for massive black holes in dwarf galaxies. Theoretical work suggests that the observed active galactic nuclei (AGNs) in dwarf galaxies are predominantly TDE-powered. To assess this claim, we perform variability analyses on the dwarf-hosted AGNs detected in the 7 Ms Chandra Deep Field-South survey, with observations spanning ≈16 yr. Based on the spectral energy distribution modeling withx-cigale, we select AGNs hosted by dwarf galaxies (stellar mass below 10¹⁰M_⊙). We focus on X-ray sources with full-band detections, leading to a sample of 78 AGNs (0.122 ≤z≤ 3.515). We fit the X-ray light curves with a canonical TDE model oft^−5/3and a constant model. If the former outperforms the latter in fitting quality for a source, we consider the source as a potential TDE. We identify five potential TDEs, constituting a small fraction of our sample. Using true- and false-positive rates obtained from fitting models to simulated light curves, we perform Bayesian analysis to obtain the posterior of the TDE fraction for our sample. The posterior peaks close to zero (2.56%), and we obtain a 2σupper limit of 9.80%. Therefore, our result indicates that the observed AGNs in dwarf galaxies are not predominantly powered by TDEs. 

ABSTRACT Broad absorption line (BAL) quasars are often considered X-ray weak relative to their optical/UV luminosity, whether intrinsically (i.e. the coronal emission is fainter) or due to large column densities of absorbing material. The SDSS-V is providing optical spectroscopy for samples of quasar candidates identified by eROSITA as well as Chandra, XMM, or Swift, making the resulting data sets ideal for characterizing the BAL quasar population within an X-ray selected sample. We use the Balnicity Index (BI) to identify the BAL quasars based on absorption of the C iv$$\lambda \, 1549$$ emission line in the optical spectra, finding 143 BAL quasars in our sample of 2317 X-ray selected quasars within $$1.5\le z \le 3.5$$. This observed BAL fraction of $$\approx$$ 6 per cent is comparable to that found in optically selected samples. We also identify absorption systems via the Absorption Index (AI) which includes mini-BALs and NALs, finding 954 quasars with AI $>0$. We consider the C iv emission space (equivalent width versus blueshift) to study the BAL outflows within the context of the radiatively driven accretion disc–wind model. X-ray selection excludes the highest outflow velocities in emission but includes the full range of absorption velocities which we suggest is consistent with the BAL gas being located further from the X-ray corona than the emitting gas. We observe both X-ray weak and X-ray strong BALs (via the optical-to-X-ray spectral slope, $$\alpha _\text{ox}$$) and detect little evidence for differing column densities between the BAL and non-BAL quasars, suggesting the BALs and non-BALs have the same shielding gas and intrinsic X-ray emission. 

Abstract Multiyear observations from the Sloan Digital Sky Survey (SDSS) Reverberation Mapping ​(RM) project have significantly increased the number of quasars with reliable RM lag measurements. We statistically analyze target properties, light-curve characteristics, and survey design choices to identify factors crucial for successful and efficient RM surveys. Analyzing 172 high-confidence (“gold”) lag measurements from SDSS-RM for the Hβ, Mgii, and Civemission lines, we find that the Durbin–Watson statistic (a statistical test for residual correlation) is the most significant predictor of light curves suitable for lag detection. The variability signal-to-noise ratio and emission-line placement on the detector also correlate with successful lag measurements. We further investigate the impact of the observing cadence on the survey design by analyzing the effect of reducing observations in the first year of SDSS-RM. Our results demonstrate that a modest reduction in the observing cadence to ∼1.5 weeks between observations can retain approximately 90% of the lag measurements compared to twice-weekly observations in the initial year. Provided similar and uniform sampling in subsequent years, this adjustment has a minimal effect on the overall recovery of lags across all emission lines. These results provide valuable inputs for optimizing future RM surveys. 

The black hole occupation fraction (focc) defines the fraction of galaxies that harbor central massive black holes (MBHs), irrespective of their accretion activity level. While it is widely accepted that focc is nearly 100% in local massive galaxies with stellar masses M⋆ ≳ 1010 M⊙, it is not yet clear whether MBHs are ubiquitous in less-massive galaxies. In this work, we present new constraints on focc based on over 20 yr of Chandra imaging data for 1606 galaxies within 50 Mpc. We employ a Bayesian model to simultaneously constrain focc and the specific accretion-rate distribution function, p(λ), where the specific accretion rate is defined as λ = LX/M⋆, where LX is the MBH accretion luminosity in the 2─10 keV range. Notably, we find that p(λ) peaks around 1028ergs−1M⊙−1 ; above this value, p(λ) decreases with increasing λ, following a power law that smoothly connects with the probability distribution of bona fide active galactic nuclei. We also find that the occupation fraction decreases dramatically with decreasing M⋆: in high-mass galaxies (M⋆ ≍ 1011−12 M⊙), the occupation fraction is >93% (a 2σ lower limit), and then declines to 66%−7%+8% (1σ errors) between M⋆ ≍ 109−10 M⊙, and to 33%−9%+13% in the dwarf galaxy regime between M⋆ ≍ 108−9 M⊙. Our results have significant implications for the normalization of the MBH mass function over the mass range most relevant for tidal disruption events, extreme mass ratio inspirals, and MBH merger rates that upcoming facilities are poised to explore.